MAXIMIZING NEUMANN FUNDAMENTAL TONES OF 

TRIANGLES 



R. S. LAUGESEN AND B. A. SIUDEJA 

Abstract. We prove sharp isoperimetric inequalities for Neumann eig- 
envalues of the Laplacian on triangular domains. 

The first nonzero Neumann eigenvalue is shown to be maximal for 
the equilateral triangle among all triangles of given perimeter, and hence 
among all triangles of given area. Similar results are proved for the 
harmonic and arithmetic means of the first two nonzero eigenvalues. 



1. Introduction 

Eigenvalues of the Laplacian arise in physical models of wave motion, dif- 
fusion (such as heat flow) and quantum mechanics, namely as frequencies, 
rates of decay and energy levels. The eigenvalues are constrained by geomet- 
ric considerations. For example, writing fj,% for the first nonzero eigenvalue 
of the Laplacian under Neumann boundary conditions on a domain of area 
A in the plane, one has that 

[i\A is maximal for disks. 

This result is due to Szego [36] for simply connected domains, with the 



extension to all domains and all dimensions by Weinberger |37j . 

Thus a free membrane of given area has highest fundamental tone when 
the membrane is circular, and the temperature of an insulated region of given 
volume will relax most quickly to equilibrium when the region is spherical. 

We prove a sharper result for triangular domains in the plane: 

\jl\A is maximal when the triangle is equilateral. 

It remains open to extend this result to n-gons (n > 4), and to find a higher 
dimensional result involving tetrahedra. 

Our result on fi±A for triangles generalizes in three different ways: to a 
stronger geometric functional, to a stronger eigenvalue functional, and to a 
trade-off between the two. To strengthen the geometric functional, we write 
L for the perimeter and prove 

H\L 2 is maximal for the equilateral triangle, 

which implies the result for \i\A by invoking the triangular isoperimetric 
inequality. 
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Strengthening instead the eigenvalue functional, we show 

/ 1 1 Y- 1 

( 1 I A is maximal in the equilateral case. 

That is, our result on the fundamental tone extends to the harmonic mean 
of the first two nonzero eigenvalues. 

We trade off a further strengthening of the eigenvalue functional against a 
weakening of the geometric functional. Specifically, we show the arithmetic 
mean (hi + ^2)/^ of the first two non-zero eigenvalues is maximal for the 
equilateral triangle, after normalizing the ratio of the square of the area to 
the sum of the squares of the side lengths. 

Our primary method is Rayleigh's Principle and the method of trial func- 
tions. Linearly transplanted eigenfunctions of the equilateral triangle are 
used to handle triangles that are close to equilateral, and linear or quadratic 
trial functions handle all the others. Neither the conformal mapping ap- 
proach of Szego nor the "radial extension" method of Weinberger seems to 
work for triangles. 

Our triangle results suggest new open problems for general domains, such 
as a possible strengthening of the Szego-Weinberger bound by an isoperi- 
metric excess term, as explained in Section [9l 

Our companion paper [22] minimizes Hi among triangles, under a diam- 
eter normalization, with the minimizer being the degenerate acute isosceles 
triangle. We know of no other papers in the literature that study sharp 
isoperimetric type inequalities for Neumann eigenvalues of triangles. The 
Neumann eigenfunctions of triangles were investigated for the "hot spots" 
conjecture, by Bahuelos and Burdzy [7j. 

Dirichlet eigenvalues of triangles have received considerable attention [TJ 
El H31 HH1 [Ml EU GSj , as discussed in Section [TOl Dirichlet eigenvalues of 
degenerate domains have also been investigated lately [HUH]- 

For a modern perspective on the Szego-Weinberger result, including its 
role as a prototype for Payne-Polya-Weinberger type inequalities, see the 
survey paper by Ashbaugh [3]. Generalizations of the Szego-Weinberger 
result to closed surfaces such as the Klein bottle, the sphere, genus 2 surface, 
projective plane and equilateral torus are known too [T2l [191 [20] [23J, [27] . For 
broad surveys of isoperimetric eigenvalue inequalities, one can consult the 
monographs of Bandle [6] , Henrot [H] , Kesavan [21] and Polya-Szego [32] . 



2. Notation 

The Neumann eigenfunctions of the Laplacian on a bounded plane do- 
main Q with Lipschitz boundary satisfy —Au = fj,u with natural boundary 
condition du/dn = 0. The eigenvalues Hj are nonnegative, with 

= no < Hi < fi 2 < ■ ■ ■ — > 00. 
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Call jUi the fundamental tone, since yT^I is proportional to the lowest 
frequency of vibration of a free membrane over the domain. Call the eigen- 
function u\ a fundamental mode. The Rayleigh Principle says 

/il = min R[v] 
u-Ll 

where w 1 1 means J n v dA = 0, and where the Rayleigh quotient is 

For triangular domains, we write: 

h > h > h > for the lengths of the sides, 
L = l\ + I2 + I3 = perimeter, 

S 2 = l\ + I 2 + Zf = sum of squares of side lengths. 
We will not need this next fact, but it is interesting that S 2 = 361 /A where 
/ is the moment of inertia of the triangular region |15l formula (6)] . 
Given nonnegative numbers a,-, define their 

aH h a n 



arithmetic mean = M{a\, . . . , a n ) 



n 



Denote the first positive roots of the Bessel functions Jo, Ji,J[, by 
jo,i * 2.4048, ~ 3.8317, ^ 1-8412. 

3. Isoperimetric upper bounds on the fundamental tone 

The Szego-Weinberger result says that among all domains of given vol- 
ume, the first nonzero Neumann eigenvalue of the Laplacian is maximized 
by a ball. Thus in two dimensions, 

ViA < Trijl,) 2 . (3.1) 

Our first theorem proves a stronger inequality for triangles. 

Theorem 3.1. For all triangles, 

ViS 2 < — (3.2) 

and hence 

lixL 2 < 16vr 2 (3.3) 

and 

Air 2 

m A<^. (3.4) 
In each inequality, equality holds if and only if the triangle is equilateral. 
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Inequality (|3.4[) for triangles improves significantly on Szego and Wein- 
berger's estimate f)3. If) . because 47r 2 /3v / 3 ~ 7.6 is much less than 7r(j'( 1 ) 2 ~ 
10.7. 

The implications (|3.2p =>- (|3.3p =>- (|3.4p are immediate from the following 
geometric inequalities. 

Lemma 3.2. For all triangles, 

UVZA <L 2 < 3S 2 . 

In each inequality, equality holds if and only if the triangle is equilateral. 

The left hand inequality is simply the triangular isoperimetric inequality. 
It implies that the triangular isoperimetric excess 

Ct = ^-A (3.5) 

is nonnegative and equals only for equilateral triangles. 
Theorem 3.3. For all triangles, 

/ IT 2 \ 4tt 2 



ill ) 3^/3 

with equality only for equilateral triangles. Equality also holds asymptotically 
for degenerate obtuse isosceles triangles. 

The discussion in Section [9] motivates such bounds involving the isoperi- 
metric excess, and shows that Theorem 13.31 implies Theorem 13.11 

So far we have maximized the fundamental tone under normalizations 
of the area, perimeter and sum of squares of the side lengths. If instead 
one normalizes the longest side, which equals the diameter of the triangle, 
then the optimal result is known already: for all convex plane domains of 
diameter D, 

ViD 2 < 4jo 2 i 

by work of Cheng [11} Theorem 2.1]. This estimate saturates for degenerate 
obtuse isosceles triangles, as discussed for example in our companion paper 
[22"1 Proposition 3.6]. 

One can bound the harmonic mean of the first two nonzero eigenvalues. 
For simply connected domains in two dimensions, the optimal inequality 
under area normalization is 

ff(Ml,M2)A<7r$,i) 2 (3.6) 

with equality for disks, by work of Szego and Weinberger [37, p. 634]. (For 
non-simply connected domains, the best result to date is H(ni, fj,2)A < 4tt 
by Ashbaugh and Benguria [4].) For triangles we have a stronger result: 

Theorem 3.4. For all triangles, 

4vr 2 

ff(/il,/X2)A<^= 
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with equality if and only if the triangle is equilateral. 

An even stronger inequality is conjectured in Section [91 using perimeter. 

Obviously Theorem 13.41 for the harmonic mean implies inequality (|3.4p 
for the first eigenvalue. 

Next, we strengthen the harmonic mean of the eigenvalues to the arith- 
metic mean, at the cost of weakening the geometric functional from A to 
A 2 jS 2 . 

Theorem 3.5. For all triangles, 

A 2 7T 2 
M(/ii,/i 2 )^2 ^ -g 

with equality if and only if the triangle is equilateral. 

By multiplying the inequalities in Theorems 13.41 and 13.51 we obtain an 
estimate on the geometric mean of the first two nonzero eigenvalues. 

Corollary 3.6. For all triangles, 

A 3 4tt 4 

- 

with equality if and only if the triangle is equilateral. 

A stronger inequality is conjectured in Section [9l using just the area. 

4. Eigenfunctions of the equilateral triangle 

This section gathers together the first three Neumann eigenfunctions and 
eigenvalues of the equilateral triangle, which we use later to construct trial 
functions for close-to-equilateral triangles. 

4.1. The equilateral triangle. The modes and frequencies of the equilat- 
eral triangle were derived two centuries ago by Lame, albeit without a proof 
of completeness. We present the first few modes below. For proofs, see the 
recent exposition (including completeness) by McCartin [26J, building on 
work of Prager [33] . A different approach is due to Pinsky [30] . 

Consider the the equilateral triangle E with vertices at (0,0), (1,0) and 
(1/2, v3/2). Then /io = 0, with eigenfunction uq = 1, and 

16^ 2 

Ml = M2 = ~9~ 

with eigenfunctions 

ui(x,y) = 2^cos [-{2x - 1)J + cos {—^)\ sm {^i 2x ~ 

U2{x,y) = cos [-^^x - 1)) - 2cos (-(2x - 1)) cos (y^J- 

Clearly u\ is antisymmetric with respect to the line of symmetry {x = 1/2} 
of the equilateral triangle, since u±(l — x,y) = —u\(x,y), whereas U2 is 
symmetric with respect to that line. 
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It is easy to check that equality holds for the equilateral triangle in The- 
orems Em 13.31 13.41 and 13.51 because m = 16tt 2 /9 and S 2 = 3, L = 3 and 
A = v/3/4. 

We evaluate some integrals of u\ and U2, for later use: 



ufdA = / u 2 , dA 
Je 



3^3 



f ( S -^pA = I (^) 2 d.4 = 32 ' r2 +- 243 
J E \ dx J J E ^dy J 32^3 

\dy J J E \ dx J 32^3 



<9x dy Je dx dy 

and also some integrals of cross-terms: 

f du\ duo , , f du\ duo , . 
u lU2 dA= / — I— icL4 = / — -— -dA = 
J E dx dx J E dy dy 

du\ du2 , . 81V3 

Ti f; — dA = h 7T 

# ax ay 32 

aui (9ti2 , . 81v3 
— — dA = 7T 

E dy dx 32 

4.2. Transplanting the eigenfunctions. Here we transplant functions 
from the equilateral triangle E to an arbitrary triangle T. Assume T has 
vertices at (—1,0), (1,0) and (a, b), where b > 0. Write 

q = a 2 + b 2 + 3. 

Let r be the affine transformation of E to T that maps the vertices (0,0), 
(l,0),(l/2,\/3/2) to (-1,0), (1,0), (a, 6), respectively. Its inverse is 

r-\x, y) = ((l+x- ay/b)/2, V3y/2b) . 

Given a function u on E, define v = u o t~ 1 on T by 

v(x, y) = (u o r _1 )(x, y) = u((l + x - ay/b)/2, V3y/2b) . 

If u has mean value zero, J B u dA = 0, then so does v, with j T v dA = 0. By 
straightforward changes of variable, 

J T \Vv\ 2 dA _ J E [(a 2 + b 2 )u 2 x - 2y/3au x u y + 3u 2 ] dA 



J T v 2 dA 4b 2 J E u 2 dA 



(4.1) 



In particular, taking a linear combination u = r yu\+5u2 of the eigenfunctions 
on E, we let v = u o r _1 to deduce 

f T \Vv\ 2 dA _ [(32vr 2 + 243)<? - 1458]7 2 - 972^a-/5 + [(32vr 2 - 243)g + 1458]5 2 
Jj, v 2 dA ~ i44& 2 ( 7 2 + 5 2 ) 

(4.2) 
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by substituting u = 7U1 + 5u 2 into (|4.ip and recalling the integrals in Sec- 
tion HU 

Similarly, putting v\ = u% o r _1 and t>2 = U2 o r _1 implies 

J T Vi?i • Vv 2 dA 486^3a 
$ T \V Vl \ 2 dA ~ "(32tt 2 + 243)g- 1458' ( ^ 

by changing variable back to E" and then using integrals from Section 14.11 

5. Proof of Theorem 13.11 and Lemma 13.21 

First we prove Lemma 13.21 Recall that M(a±, . . . ,a n ) denotes the arith- 
metic mean and define 

G(a%, . . . , a n ) = tya\ ■ ■ ■ a n = geometric mean, 
Q(ai, . . . , a n ) = \J (af + • • • + a%)/n = quadratic mean. 
Then for any triangle, 

12y/3A = 3V3LG(L - 2/i, L - 2l 2 ,L- 2/ 3 ) 3/2 by Heron's formula 

< 3v / 3LM(L — 2l\,L— 2/ 2 , L - 2hf' 2 

= L 2 

= 9M(h,l 2 ,l 3 ) 2 

<9Q(h,l2,kf 
= 3S 2 . 

Inequalities between these means become equalities if and only if l\ = l 2 = h, 
meaning the triangle is equilateral. Thus Lemma 13.21 is proved. 

For Theorem 13. H it remains to prove H\S 2 < 167r 2 /3 with equality if and 
only if the triangle is equilateral. 

Let T be a triangle. By rescaling and rotating and reflecting, we can 
assume the longest side of T has length 2 with vertices at (—1,0) and (1,0), 
and that the third vertex (a, b) satisfies a > and b > and 

(a + l) 2 + 6 2 <2 2 . (5.1) 

Let us express these constraints in terms of new variables 

p = b 2 and q = 3 + a 2 + b 2 . 

By definition, q > 3 and < p < q — 3. Condition (|5.ip says 2a < 6 — q, and 
since a > we conclude q < 6. Hence 

3<g<6, 0<p<q-3. (5.2) 

Further, by substituting a = y/q — 3 — p into 2a < 6 — q and then squaring, 
we find 

p>Q c (q) = ~q 2 +4q- 12. (5.3) 
The constraint region determined by (|5.2p - (|5.3p is plotted in Figure [TJ 




Figure 1: The constraint region in the proof of Theorem 13, 1\ showing the 
line p = q — 3 and the curve p = Q c (q)- 



When q = 6 the constraints require p = 3, so that a = 0, b = y/3, 
and so T is equilateral. In that case equality holds in the theorem, with 
HiS 2 = 167r 2 /3 by Section [U So from now on we assume q < 6. 

Continuing with the proof, note the squares of the side lengths of the 
triangle add up to 

S 2 = 2 2 + (a - l) 2 + b 2 + {a + l) 2 + b 2 = 2q. (5.4) 

Consider now the linear functions 

f{x, y)=x-^, g(x, y)=y--, 

which integrate to zero over the triangle T. Our first trial function is the 
linear combination v = f + 75 where 7 G R. By Rayleigh's Principle, 

MlS 2 <fl[/ + 79]S 2 

= ^L±^ 2 „ (5 . 5) 
3 + (a + jb) z 

This last expression is less than 167r 2 /3 (as desired for the theorem) if 

27(1 + 7 2 )g - 4vr 2 (3 + (a + 7 &) 2 ) < 0. 

The left hand side is a quadratic polynomial in 7, and hence an appropriate 
7 exists if the discriminant is positive, which is equivalent to 

P<Qi(q) = -^(^ 2 -27)q 2 . (5.6) 

For our second trial function, let u\ be the first antisymmetric eigenfunc- 
tion of the equilateral triangle E and recall the affine transformation r from 
E onto T, as described in Section HI The transplanted function v = u\ or -1 
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integrates to over T, and so can be used as a trial function. By Rayleigh's 
Principle, 

^S 2 < R[ Ul or-^S 2 

(32?r 2 + 243)q - 1458 

where the Rayleigh quotient has been evaluated by formula (|4.2p with 7 = 1 
and 5 = 0. Notice the last expression is less than lQir 2 /3 if 

1 

384^r 2 

Observe Q c (q) > Qi(o) when 



P > Qi(<?) = rr-^((32vr 2 + 243) 9 2 - 1458g). (5.7) 



768vr 2 

5.03 ~ = < q < 6. 

128vr 2 + 243 H 

Therefore if 5.04 < q < 6 then the constraint (|5.3j) implies p > Q c (q) > 
Qi(q), so that (|5.7p holds and hence //iS 12 < 167r 2 /3. 
Next, observe Qt(q) > Q\{q) when 

1458tt 2 

° <9< 32vr4 -621^ + 5832 " 5 - 10 - 
Therefore if 3 < q < 5.09 then for all p E R, either p < Qe(q) or else 
P > Q^(g) > Qi(q)i so that either f)5.6j) or (|5.7p holds; in either case, we 
conclude ^iS* 2 < 167r 2 /3. 

We have proved /UiS 2 < 167r 2 /3 in the whole constraint region 3 < q < 6, 
and so the proof is complete. 

To summarize the above proof, notice that close-to-equilateral triangles 
(with 5.04 < q < 6 above) are handled by the trial function u\, which is 
the linearly transplanted eigenfunction of the equilateral triangle, while far- 
from-equilateral triangles (3 < q < 5.09) are treated with either that same 
transplanted eigenfunction or else the linear trial function / + 75. 

6. Proof of Theorem [33] 

Equality in the theorem holds for equilateral triangles, as observed in 
Section^ In addition, when an obtuse isosceles triangle degenerates towards 
a line segment, equality holds in the limit because A — > 0, L — > 2D and 
w _> 4j 2 1 L>" 2 (by [221 Proposition 3.6]). 

Assume for the rest of the proof that the triangle is non-equilateral. By 
rescaling, rotating and reflecting, we reduce to considering the triangle T 
with vertices (—1,0), (1,0) and (a, 6), where a > 0, b > and all the side- 
lengths are less than or equal to 2. This triangle has area A = b and diameter 
D = 2. 

Introduce new parameters r = (I2 + h)/2 and s = (I2 — fa)/2 defined in 
terms of the sidelengths 



l 2 = v^a + l^+ft 2 and l 3 = ^(a-l) 2 + b 2 . 
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These new parameters occupy a triangular region in the rs-plane (see Fig- 
ure [2] below) with 

Kr < 2, < s< 1, r + s < 2, 

since ?3 < ?2 < ^1 = 2 and £1 < I2 + ^3 and /2 < ^1 + ^3- The line r = 1 
corresponds to degenerate triangles. Since T is not equilateral, we know 
(r,s) ^(2,0). 

In terms of the new parameters, we have 



a = 


rs, 




6 2 = 


(r 2 -l)(l- 




9 = 


a 2 + 6 2 + 3 


= r 2 + s 2 + 2, 


L = 


^1 + k + ^3 


= 2(l+r), 


A = 


b = V(r 2 - 


-1)(1- S 2 ), 


£t = 


T 2 

A 

12^/3 


3^/3 



-l)(l-,2). 



To handle close-to-degenerate triangles, we recall Cheng's bound H\D 2 < 
^Jqi for convex domains (see Theorem 2.1]); since our triangle T has 
diameter D = 2, Cheng's bound gives /xi < j 2 1 and thus 

2 

Ml (>1 + ^-^t) < JliA + tt 2 £ T - 
v Jo,i J 

This last expression is less than 47r 2 /3-\/3 (as desired for the theorem) if 

(r - l)( r + 3) - (l - V27(r 2 - 1)(1 - s 2 ) < 0. (6.1) 

The left side is increasing with s. Putting s = 2 — r (the largest value of s 
in our parameter region), we find that (|6,ip holds if 

(r + 3) - (l - ^) v / 27(r + l)(3-r) < 0. 

The expression on the left is convex for — 1 < r < 3, and is negative at r = 1 
and r = 5/4, and so inequality (|6.1j) certainly holds for 1 < r < 5/4. Thus 
the theorem is proved in that range, as indicated in Figure [2j Next, suppose 
5/4 < r < 3/2 and < s < 1/3. Replacing s by 1/3 in (|6.ip . and replacing 
r + 3 by (3/2) + 3 = 9/2, we see that it suffices to prove 

(r " I)©' " (l " fr)' • 27 (r + 1)(1 - (1/3) 2 ) < 0. 

This linear inequality is easily established when r < 3/2. Hence the theorem 
holds when 5/4 < r < 3/2 and < s < 1/3, as indicated in Figure [2j 
For the remaining parameter regions in Figure [21 we will show 

R M (A + f (6.2) 
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s * 




Figure 2: Different shadings represent different trial functions used to prove 
Theorem 13.31 The darkest shading represents the circular sector-based trial 
function underlying Cheng's bound; intermediate shading represents a linear 
trial function; and the lightest shading represents a trial function based on 
eigenfunctions of the equilateral triangle. The point (2, 0) corresponds to 
the equilateral triangle, and degenerate triangles have r = 1. 



for some trial function v having mean value zero; notice here we have re- 
placed 7T 2 /_7q X with the slightly larger value 12/7. 

First, consider the linear trial function v(x,y) = x — a/3, which has 
mean value zero and R[v] = 18/(3 + a 2 ). The desired estimate (16. 2h is then 
equivalent to 

(r + l) 2 - ^(3 + rV) < A^ 27(r 2_ 1)(1 _ s2) . (6 . 3) 

Restrict to the region where 5/4 < r < 3/2 and 1/3 < s < 2 — r. The left 
and right sides of inequality ()6.3p each decrease as s increases, and so we 
put s = 2 — r in the right side and s = 1/3 in the left, reducing our task to 
proving the inequality 

(r + I) 2 " ^(27 + r 2 ) - A V27(r 2 - l)(r - 1)(3 - r) < 0. 

The expression is convex for 1 < r < 3, and so it is enough to verify the 
inequality at the endpoints r = 5/4 and r = 3/2. Direct calculation shows it 
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is true at those endpoints, and so the theorem is proved when 5/4 < r < 3/2 
and 1/3 < s < 2 - r. 

Next restrict to the region 3/2 < r < 8/5 and 2/5 < s < 2 — r. Like above, 
we put s = 2 — r in the right side and s = 2/5 in the left. The resulting 
expression is again convex and the inequality is true at the endpoints r = 3/2 
and r = 8/5. 

The third case uses a stronger version of (16. 2ft that is easier to handle. 
We have 

. 12 L 2 5 L 2 180A 2 

A+ T £T= W = 3~~7 A ^W = 3~^T 2 
by the triangular isoperimetric inequality 12\^3A < L 2 . Thus (|6.2p will hold 
if we can show 

A 2 287T 2 

R[v}(L 2 _180 IJ )<— . (6.4) 

Consider the eigenfunctions ui and «2 of the equilateral triangle -E, as 
in Section U and the affine transformation r from E to T. Transplant the 
eigenfunctions by v\ = uiot~ 1 and v 2 = U2°t~ 1 , so that v\ and V2 integrate 
to over T. Then take the trial function 

1 

V = Vl + -v 2 . 

Its Rayleigh quotient R[v] can be evaluated by formula f|4.2j) with 7 = 1 and 
5 = 1/3. Substituting in this formula for R[v] reduces the desired estimate 

dS3D to 

U(r, s)V (r, s) + W(r, s) < (6.5) 

where 

Ufa s ) = f^l + 216j (r 2 + s 2 ) - 324 v / 3rs + - 864, 

Vfa s) = 4(r + l) 3 + 45(r - l)(s 2 - 1), 

Wfa s) = 448Q7r2 (r - l)(r + 1) V - 1). 
3 

We will prove inequality (|6.5p in the trapezoidal region defined by 3/2 < 
r < 2 and < s < min {2/5, 2 — r}. First we show UV + W is convex with 
respect to s. The second derivative of UV + W with respect to s is 

9 o N 6407T 2 , 9 , N 

432(32 - 132x + 114x 2 + 49x 3 ) + — - — (32 + 306a; + 282x 2 + 91x 3 ) 

9 

-87480a;(x + l)y/3s + 480(40vr 2 + 243)xs 2 , 

where we have put x = r — 1 > 0. We want to show this quantity is positive. 
We may discard the s 2 term, since it is certainly positive. Further, we may 
replace \^3s by the larger number 1, noting s < 2/5 in the trapezoidal region, 
and we may replace ir 2 by the smaller number 9. After these reductions, 
one is left with the polynomial 

34304 + 51336x + 142248x 2 + 79408x 3 , 
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which is positive. Hence UV + W is convex with respect to s, in the trape- 
zoidal region. 

By convexity, it suffices to prove (16. 5j) on the upper and lower boundary 
portions of the trapezoid, that is, where 3/2 < r < 2 and s = 0, or 3/2 < 
r < 8/5 and s = 2/5, or 8/5 < r < 2 and s = 2 — r. 

We start with 3/2 < r < 2 and s = 0. By substituting s = and 
r = x + 3/2 into (|6.5p . we reduce to the inequality 

- l/2)((160vr 2 + 972)x 4 + (1760tt 2 + 10692)x 3 + (4680vr 2 + 32805)x 2 

+(3400vr 2 + 35235)x + (-3100tt 2 + 34020)) < 0, 

which obviously holds true for < x < 1/2, that is, for 3/2 < r < 2. 
Incidentally, the root at x = 1/2 arises from the equilateral triangle (r = 2). 

Now we take 8/5 < r < 2 and s = 2 — r. By substituting s = 2 — r and 
r = x + 8/5 into (|6.5p . we reduce to the inequality 

(x - 2/5)((7000000vr 2 + 7441875c)x 4 + (21400000tt 2 + 15278625c)x 3 
+(34200000tt 2 + 8693325c)x 2 + (7976000tt 2 + 12510855c)x 

+ (-10211200vr 2 + 11572632c)) < 0, 

where c = 4 + 3\/3. Note the constant term in the quartic is positive. Thus 
the inequality holds when < x < 2/5, that is, when 8/5 < r < 2. 

Lastly we take 3/2 < r < 8/5 and s = 2/5. By substituting s = 2/5 and 
r = x + 3/2 into (|6.5p . we reduce to the inequality 

(243000 + 40000?r 2 )x 5 + (2551500 - 145800\/3 + 420000vr 2 )x 4 
+(7341030 - 1312200\/3 + 1095600vr 2 )x 3 + (6530625 - 2996190\/3 
+934600vr 2 )x 2 + (4352859 - 3623130\/3 - 138480vr 2 )x 

+(-4211433 - 2383830^3 + 820260^) < 0. 

The coefficients of the second and higher powers of x are positive, and so 
this quintic polynomial is convex. The polynomial is negative at x = and 
x = 1/10, and hence is negative whenever < x < 1/10, that is, whenever 
3/2 < r < 8/5. This observation completes the proof. 

7. Proof of Theorem S3] 

The harmonic mean of the first two nonzero eigenvalues is characterized in 
terms of the Rayleigh quotient by Poincare's variational principle [BJ p. 99]: 

H{^^ 2 ) = rffin {#(#[/!], i?[/ 2 ]) : f x J_ 1, / 2 J_ 1 and V/i + V/ 2 }. 

As in the earlier proofs, we reduce to considering the triangle T with 
vertices (-1,0), (1,0) and (a, b) where a > 0, b > and (a + l) 2 + b 2 < 2 2 . 
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Notice a < 1 and b < \/3, and that the triangle has area A = b. Recall the 
definition 

q = a 2 + b 2 + 3 

and observe 

q > b 2 + 3 > 2V3b. (7.1) 
Consider the polynomial trial functions 



fi(x,y) = x -~, 
X ~ 3/ 



3 + a 2 
18 



whose coefficients have been chosen to ensure /i and f 2 have mean value 
zero over T and have orthogonal gradients (J* T V/i • Vf 2 dA = 0). These 
trial functions have Rayleigh quotients 

18 „ r . , 360 



3 + a 2 ' UZJ 7(3 + a 2 ) 

Hence by Poincare's principle, 

H{iii,iJi2)A<H{R[h},R[f 2 ])A 
80b 



3(3 + a 2 ) 

This last expression is less than Air 2 / 3^/3 if 



b+^f<l- (7.2) 



20^/3 £ 
T 2 " < 6' 

Thus strict inequality holds in the theorem if (|7.2[) is true. 

Next consider the eigenfunctions u\ and u 2 of the equilateral triangle E, 
as in Section 01 and the affine transformation r from E to T. Transplant 
the eigenfunctions by v\ = u\ o t _1 and V2 = «2 t _1 , so that and ^2 
integrate to over T. Let 7 € R. The trial functions V\ and 7U1 + ^2 have 
Rayleigh quotients 

(32tt 2 + 243)<? - 1458 



1446 2 

[(32tt 2 + 243)g - 1458]7 2 - 972^7 + [(32tt 2 - 243)g + 1458] 

1446 2 ( 7 2 + 1) 



R[vi] 
R[yvi + v 2 ] 

as shown by formula (|4.2p with 5 = and <5 = 1, respectively. 

We choose the coefficient 7 so that the gradients of v± and 7^1 + v 2 are 
orthogonal: 

L Vui • Vv 2 dA 486 y/3a 

7 ~ ' " 



J T \Vvi\ 2 dA (32^ 2 + 243)g - 1458 
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by formula (14. 3p . Then by Poincare's principle, 

H(ji lif Mi)A < H(R[vi],R[yvi + v 2 ])A 

_ (1024vr 4 - 243 2 )g 2 + 12 • 243 2 6 2 
~ 4608vr 2 6g ' 

which is less than 47r 2 /3\/3 if and only if 

r q 2^3 • 243 2 , , 

Hence when (\7.3\i is true, strict inequality holds in the theorem. 

To complete the proof, we divide into three cases. First, if q/b = 2^3 
then T is equilateral (because a = and 6 = y/3 by considering equality in 
(|7.ip ). so that equality holds in the theorem by Section [H Second, if 

2^3 < 7 < 5 
o 

then strict inequality holds in the theorem by (|T.3j) . Third, suppose 

which means 56 < a 2 + b 2 + 3. Since a < 1 we deduce < (b — 1) (6 — 4), so 
that b < 1. Therefore estimate (|7.2p is true, because its left side is at most 
1 + 20\/3/V 2 ~ 4.51 while its right side is at least 5. Hence once again the 
theorem holds with strict inequality. 

8. Proof of Theorem GTS 

The arithmetic mean of the first two nonzero eigenvalues is character- 
ized in terms of the Rayleigh quotient by Poincare's variational principle [H 
p. 98]: 

M( m ,/i 2 ) = mm{M(fl[/i],l?[/2]) : /i J_ l,/ 2 J_ 1 and f x _L / 2 }. 

Like in the earlier proofs, we need only consider the triangle T with ver- 
tices (—1,0), (1,0) and (a, 6), where a > 0,6 > and 

(a + l) 2 + 6 2 < 2 2 . (8.1) 

Recall the eigenfunctions u\ and n 2 of the equilateral triangle E, as in Sec- 
tion HI and the affine transformation r from E to T. Transplant the eigen- 
functions by vi = u\ o t -1 and v 2 = u 2 o r _1 , so that the trial functions t?i 
and V2 integrate to over T. Note j T v\V2 dA = by the antisymmetry and 
symmetry properties of u\ and u 2 . 
The Rayleigh quotients evaluate to 

„. . (32vr 2 + 243W - 1458 , „. . (32vr 2 - 243)g + 1458 
R[V1] = 144P ^ R[V2] = 144P 
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by formula (14. 2p with 8 = and 7 = 0, respectively. Hence 

M(/ii,/i 2 ) <M(R[ Vl ],R[v 2 ]) 
2ir 2 q 
9b 2 
ir 2 S 2 



9A 2 

since S 2 = 2q by formula (|5.4p . and A = | • 2 • b = b. This last estimate is 
the desired inequality. 

If the triangle is equilateral then equality holds in the theorem, by Sec- 
tion m 

Suppose equality holds in the theorem. We will show T is equilateral. 
Since equality holds in our argument above, the arithmetic mean of and 
fj-2 equals M(R[vi], R[v2]), which implies by the proof of the variational 
principle (see [61 p. 98]) that the span of v\ and v 2 equals the span of some 
two eigenfunctions with eigenvalues \x\ and [i% these eigenfunctions can 
be assumed orthogonal in L 2 (T). Hence there exists a linear combination 
v = jvi + Sv2 (with coefficients 7, 8, not both zero) that is an eigenfunction 
of the Laplacian on T. By direct calculation, 

Av 2ir 2 , 9 , 9 N , 

— = 7^ a +b 9 at x, y = 0,0, 
v 9b z 

Av 2ir 2 2 2 

— = -— y(a +6+3) at (x,y) = (a,b), 
v 9b z 

where we used the definition v = (7^1 + 8112) or -1 and called on the formulas 
for m, U2 and r _1 in Section HI These expressions for (Av)/v at (0, 0) and at 
(a, 6) must be equal, since (Av)/v is constant by the eigenfunction property. 
Hence a 2 +b 2 = 3. The constraint (|8.ip then implies a < 0, so that a = and 
hence 6 = \/3- Thus T is equilateral, completing the proof of the equality 
statement. 

9. Discussion of isoperimetric excess, and open problems 



One could attempt to strengthen Szego and Weinberger's result (|3.ip for 
general plane domains by adding a multiple of the general isoperimetric 
excess, which is defined by 

L 2 

£ = — -A. 

47T 

Note the excess is nonnegative by the isoperimetric inequality, and that it 
equals only for disks. 

Problem 9.1. Does there exist 8 > such that 

for all convex plane domains? What is the largest possible 8? 
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Among general domains, 

\x\l? is not maximal for the disk, 

because the equilateral triangle and the square both have [i\L 2 = 16tt 2 ~ 
158, which exceeds the value 4-7r 2 (j'J x ) 2 ~ 133 for the disk. Hence Prob- 
lem [9J] needs 5 < 1, because when 5 = 1 one has A + 5£ = L 2 /4ir. 

Problem 9.2. Determine the maximizers for /J-iL 2 , among all bounded con- 
vex domains in the plane. 

The convexity hypothesis eliminates domains with fractal boundary, for 
which L is infinite and [i\ can be positive [29J. A result somewhat similar to 
Problem 19.11 was proved by Nadirashvili |28j . but with a measure theoretic 
"asymmetry" correction instead of the isoperimetric excess. 

The triangular version of Problem 19.11 is to find 5 > such that 

4vr 2 

m .(A + ( $S r )<— (9- 1 ) 

where we recall the triangular isoperimetric excess £t = (L 2 /12\/3) — A 
defined in (|3.5p . 

We have already proved triangular excess bounds of the form f)9. 1 j) : the 
perimeter bound \i\L 2 < 16ir 2 in Theorem 13.11 has that form for 6 = 1, 
because A + £ T = L 2 /12\/3. Theorem ECS is even stronger, for it proves 
(|9.ip with 5 = 7t 2 /j'q i ~ 1.7 and hence with all smaller values of 5 too, such 
as 6 = 1 and 5 = 3/2. 

Theorem 13.31 implies Theorem 13. 1\ because (I9.ip with 5 = 3/2 implies the 
sum-of-squares bound fiiS 2 < 167r 2 /3 by Lemma 19.31 below. 

Lemma 9.3. For all triangles, 

,2 . 12 I . 3 



with equality for equilateral triangles and asymptotic equality for degenerate 
acute isosceles triangles. 

Proof. In the notation of Section [6] we have 

L = 2(l + r), 



a = vV-iXi-* 2 ), 

S 2 = 2(2 + r 2 + s 2 ), 

where the parameters satisfy l<r<2,0<s<l. Substituting these 
quantities into the lemma, we see the task is to prove 



2r - (1 + s 2 ) > ^/3(r 2 -1)(1- s 2 ). 
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The left side is positive. By squaring both sides and rearranging, we reduce 
to the equivalent inequality 

/ 1 + S 2 \2 o(l-S 2 ) 2 

which is clearly true. Equality holds when r = 2, s = 0, which is the 
equilateral case. Equality holds asymptotically when s = l,r = 1, which 
corresponds to a degenerate acute isosceles triangle. □ 

Incidentally, we settled on the choice of 5 in Theorem 13.31 by increasing 5 
until some non-equilateral triangle also gave equality in the theorem. Any 
further increase would prevent the equilateral triangle from being optimal. 

Is there a lower excess bound for fii, complementing the upper bounds 
in Theorem 13.31 and Problem 19. IP 

Problem 9.4. Is there a constant 5 > such that for all triangles, 

4vr 2 

M1 (A + ^ T )>— ? 
Is there a constant 5 > such that for all bounded Lipschitz plane domains, 

For triangles one would need 5 > 4-7r 2 /j 2 1 at least, in order for the in- 
equality in the Problem to hold for the degenerate acute isosceles triangle 
(see |22[ Corollary 3.5]). For general domains one would need 5 > (j[ x ) 2 at 
least, in order for the inequality to hold for the degenerate rectangle. 

Turning now from the fundamental tone to the harmonic mean of the first 
two eigenvalues, we raise: 

Conjecture 9.5. For all triangles, 

H(fi ht i 2 )L 2 < 16vr 2 

with equality if and only if the triangle is equilateral. 

This conjecture would be stronger than Theorem 13.41 where we used A 
instead of L 2 . 

For the geometric mean our numerical work similarly suggests: 
Conjecture 9.6. For all triangles, 

,2 !6vr 4 
Mi^2-4 < 

with equality if and only if the triangle is equilateral. 

This conjecture would be stronger than Corollary 13.61 where we had 
A 3 /S 2 instead of A 2 , in view of Lemma 13.21 It would also be stronger 
than Theorem 13.41 which uses the harmonic mean and A. 
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To contrast the last two conjectures, notice that when the eigenvalue 
functional is strengthened from the harmonic mean to the geometric mean, 
the scaling factor is weakened from the perimeter to the area. 

For the geometric mean on general plane domains, Iosif Polterovich has 
conjectured f^i^A 2 < n 2 (j[ i) with equality for the disk (see Mathematis- 
ches Forschungsinstitut Oberwolfach MFO Report 6/2009). This inequality 
is known up to a factor of 2, by combining the Szego- Weinberger inequality 
Hi A < n(j'i x ) 2 with the recent result of Girouard et al. that H2A < 2ir(j[ x ) 2 
(with equality holding for a domain degenerating suitably to two disjoint 
disks of equal area) [17] . Incidentally, the better bound [12A < n(j'ii) 2 
holds for domains with 4- fold rotational symmetry [JJ. 

Lastly, for curved surfaces we raise the open problem of maximizing 
the Neumann fundamental tone among spherical and hyperbolic triangles 
of given area, in the two dimensional sphere and hyperbolic disk respec- 
tively. Note the Szego-Weinberger inequality for general domains has been 
extended from euclidean space to curved surfaces in two dimensions [H 
§111.3.3], and to the sphere and hyperbolic space in all dimensions [5] IIP]. 



10. Survey of Dirichlet eigenvalue estimates 

We close the paper by mentioning analogous results for Dirichlet eigen- 
values. The Dirichlet analogue of the Szego-Weinberger bound (|3.ip is the 
Faber-Krahn inequality 

Al^4 > 7Tio,l, 

which holds with equality for the disk. The triangular version of this in- 
equality appears in the book of Polya and Szego [32, p. 158]: 

XlA> -vs 

with equality for equilateral triangles. 

Faber-Krahn type bounds are necessarily one-sided, because a long, thin 
domain can have fixed area and Ai arbitrarily large. To obtain a two-sided 
bound on Ai one must weaken the geometric functional. For example, for 
convex plane domains one has 

vr 2 x A 2 it 2 

16 ~ T 2 ~ T' 



where the upper bound is due to Polya [31] and the lower bound to Makai 
|25j . Equality holds asymptotically in these bounds for degenerate circular 
sectors and degenerate rectangles, respectively. For triangles a sharper upper 
bound was proved by Siudeja [34] : 



vr 2 A 2 ir 2 , , 

t- < A 1T¥ < — (10.1) 



16 - A1 I^ 

with equality in the upper bound for equilateral triangles. 



20 



These last bounds can be strengthened to include the isoperimetric excess; 
see Siudeja [34], Conjecture 1.2] and Freitas and Antunes [I]. 

The geometric functional A 2 / S 2 that we combined in Theorem 13.51 with 
the arithmetic mean of the Neumann eigenvalues has been studied also in 
the Dirichlet case. Freitas |13] showed for arbitrary triangles that 

x A 2 7T 2 

which is slightly weaker than fllQ. 1H : quadrilaterals have been studied too 
|16j . Conjectures involving Ai and geometric functionals have been raised 
by Antunes and Freitas pQ. 

The Dirichlet gap conjecture for triangles, due to Antunes and Freitas [2] 
Conjecture 4], claims that (X2 — Xi)D 2 is minimal for the equilateral triangle. 
Some progress has been made recently by Lu and Rowlett [24]. 

Finally, recall the inverse problem of determining the shape of a triangle 
from knowledge of its Dirichlet spectrum. The spectrum is tremendously re- 
dundant, since it is determined by merely three parameters (the side lengths 
of the triangle). It is plausible that the triangle could be determined (up 
to congruence) by knowing just finitely many eigenvalues. Chang and De- 
Turck [9] did so nonconstructively, with the required number of eigenvalues 
depending on Ai and A2- A constructive approach or explicit formula for 
solving the inverse problem would be most welcome. 
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